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We study the coherent storage and retrieval of a very short multimode light pulse in an atomic 
ensemble. We consider a quantum memory process based on the conversion of a signal pulse into 
a long-lived spin coherence via light matter interaction in an on-resonant A -type system. In order 
to study the writing and reading processes we analytically solve the partial differential equations 
describing the evolution of the field and of the atomic coherence in time as well as in space. We 
show how to optimize the process for writing as well as for reading. If the medium length is fixed, 
for each length, there is an optimal value of the pulse duration. We discuss the information capacity 
of this memory scheme and we estimate the number of transverse modes that can be stored as a 
quantum hologram. 

PACS numbers: 42.50.Gy, 42.50.Ct, 32.80.Qk, 03.67.-a 



I. INTRODUCTION 

Storage and read out of quantum states of light in mat- 
ter represent a major challenge for quantum communi- 
cations and quantum information processing. Since the 
first proposals for such process about a decade ago (T], Q , 
and the first experimental implementations 0, |J a num- 
ber of schemes have been studied d, 0] . The main objec- 
tive is to store and retrieve light pulses or photons with- 
out destroying their quantum state, that is the system 
should be able to store at the same time two on commut- 
ing variables like the two quadratures of a light pulses 
and allow for their retrieval. Mapping quantum states 
of light like single photons or qubits [^-Q , coherent light 
pulses without noise added fToj and squeezed light Tl|, 12 1 
onto long lived states of atomic coherence has been first 
experimentally investigated in alkali-metal gases. More 
recently quantum memory for entangled photons have 
been demonstrated in rare-earth doped crystals [l^, • 

While the principle of such memory registers has been 
mainly developed for a single temporal and spatial mode 
jl5l - l2g | , it appears that the need for multiplexing will be 
high and that multimode memories need to be developed. 
Several schemes for a spectral multi-mode memory have 
been proposed based on off resonant Raman interaction 
|2l| - |23| . controlled reversible inhomogeneous broadening 
(CRIB) [2I-I23 and atomic frequency comb (AFC) 
with significant recent experimental advances in vapors 
[29I - I32I and crystals [SS^I. In the direction towards spa- 
tial multimode storage of a quantum image, a quantum 
hologram scheme was proposed in reference ^J], based 
on quantum non demolition measurement (QND) type 
interaction and in reference (ssj based on Raman-type 
interaction in non-coUinear field configuration. Experi- 
mentally storage of a classical image was demonstrated 



in an atomic vapor based on electromagnetically induced 
transparency (EIT) type interaction [3(|. In this pa- 
per, we propose a different storage protocol which signif- 
icantly combines both spectral and spatial types of mul- 
tiplexing. Our scheme provides a storage of a quantum 
image carried by a very short broadband light pulse in 
an atomic medium. We show that it allows both spectral 
and spatial multimode storage with a good efficiency. 

We consider light storage in a three-level medium in 
a A-configuration, with a strong driving field close to 
resonance with one of the transitions and a weak signal 
field close to resonance with the other transition. The two 
lower levels are assumed to be long-lived, being sublevels 
of the ground state. 

The quantum signals carried by the signal pulse will 
be stored in the long lived ground state coherence of the 
atomic ensemble. Contrary to other schemes, based on 
EIT, we consider very short pulses, actually much shorter 
than the excited state decay time. Interaction with the 
medium is then very fast and it does not allow for the 
build up of EIT. In contrast to echo-type memory, we 
consider the simultaneous action of both signal and driv- 
ing fields on the medium. 

Storage based on off-resonant Raman interaction was 
recently proposed and demonstrated for a light pulse 
much broader than the natural linewidth (23l . Isoj . The 
main difference with our protocol is that we propose to 
explore a resonant interaction and rely on more efficient 
light-mater coupling. Contrary to CRIB and AFC based 
schemes our protocol does not deal with an inhomoge- 
neous broadened medium and does not use any rephasing 
process. However, we will see that the quantum informa- 
tion storage can be quite efficient, in agreement with |T6| . 

The article is organized as follows. In SecHHthe physi- 
cal model of a 3D-memory process based on the resonant 
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FIG. 1: Three level atomic system interacting with driving 
field Q, and signal field a. 



interaction of the three-level atoms with short pulses is 
discussed in detail. In Sec. IIIII the writing of a weak 
quantum field in space and in time is considered. The 
efficiency and optimization of writing are analyzed. In 
Sec.|IV]the read-out is investigated in both forward and 
backward processes. The role of the diffraction and the 
number of stored modes are discussed at the end of this 
section. In App. I^the mathematical aspects of the prob- 
lem are considered in detail. 



II. PHYSICAL MODEL 

In this paper, we will study the ability of a system 
based on ensembles of three-level atoms (Fig. [T]) to store 
temporal as well as spatial multimode quantum fields, 
thus implementing a quantum hologram. For this, we 
take the driving field as a plane wave, while the signal 
field has as a transverse structure. A similar problem 
has been considered in Ref [13] in the case of a QND 
interaction of light with atoms. 

The signal and driving pulses are assumed to be much 
shorter than the excited state lifetime 7""'^, so that we can 
neglect spontaneous emission during the writing process, 
which eliminates a source of fluctuations and dissipation. 
For the writing process, the driving and signal pulses 
are assumed to be simultaneous and to have the same 
duration Tw- Similarly the driving pulse in the process 
of read-out has a duration Tji much smaller than 7~^. 
Moreover, we assume that duration of the pulses is much 
larger than the time it takes for the light to go through 
the atomic medium of length L 



7 ^ > Tw,R > L/c. 



(2.1) 



The atom-field interaction will be written in the dipole 
approximation, and the corresponding interaction Hamil- 
tonian reads 



(2.2) 



where r} represents the spatial coordinates of the j-th 

atom and where the field E{fj) is a combination of two 
fields, signal and driving fields, which are interacting on 
the two neighboring atomic transitions 



The driving and signal fields are pulses of equal duration 
that verify Ea. (|2.ip . The driving field is in a monochro- 
matic coherent state with frequency Ud and propagates as 
a plane wave with wavevector fc^. Since most of the atoms 
are in level |1 >, as explained below, we assume that we 
can neglect the driving field absorption and that it prop- 
agates along the z-axis through the atomic medium with 
a constant amplitude E^: 



Edir,t) = Eoe 



-iuiat+ik^z 



(2.4) 



We treat the signal field in the paraxial approximation 
as quasi-monochromatic transverse multimode wave of 
frequency LOg propagating in the same direction as the 
driving field with an average wavevector . We can then 
write the signal field as: 



Es{r,t) 



huj 



2£nC 



^^ — iijJst-\-iksZ 



a{z,p,t), (2.5) 



where a(z, p, t) is the annihilation operator for the signal 
field and p = p{x, y) describes the transverse signal field 
distribution. The amplitude a{z,p,t) is normalized so 
that the mean value {a^ {z, p,t)a{z, p,t)) is the photon 
number per second per unit area. 

Using the rotating wave approximation we can write 
the Hamiltonian as 

V = j dz (fp ihga{z,p,t)azi{z,p,t)e^^^^ (2.6) 
-ihga) (z, p, i)(7i3(z, p, t)e-*^^^ 

+ inna32{z, p, Oe^^^^ - ihn*a23{z, p, t)e"*^rf^] , 

where aik are the atomic coherence operators between 
levels i and k and dik are the corresponding dipole matrix 
elements; O = Eod23/h is the Rabi frequency for the 
driving field ; g is the coupling constant between the 
signal field and atom in the dipole approximation: 
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2eohc 



1/2 



(2.7) 



E{f,t) = E,{r,t)+Edif,t). 



(2.3) 



For the sake of the simplicity we assume d^fc to be real 
so that g = g* and we have set the frequency detunings 
(lUs — Wis) and {ujd — W23) equal to zero. 

The collective atomic coherences are given by a sum of 
individual atomic operators 

d-mn{r,t) = ^ |m >< n\jd^{f- fj) 
j 

TO,n=: 1,2,3, m n, f={z,p}. (2.8) 
These operators obey the commutation relation 

['T™„(f,i),a„™(r',t)] = (7V„(r,t)-7V„(f,i)) ,53(r-r'), 

(2.9) 
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where Nm are the cohective atomic population operators and we make the following substitutions 



N^{r,t) = ^ |m X m\j5\T- f,). (2.10) 



In the paraxial approximation, the slow field amplitude 
a{f,t) obeys the commutation relation 



id 18^ 
kg dz 2k1 dp^ 



cS^ir-r'). 
(2.11) 



In the following we will assume spatially slowly varying 
atomic operators and make the substitution 



0'12 



e 



-i(kd - ks)z 



0'12- 



(2.12) 



The driving and signal fields are assumed to be super- 
imposed at all times. They start interacting with the 
atomic medium at time t — 0. The input plane of the 
atomic medium is located at z = 0. We assume that all 
the N atoms are initially in state |1). 

The complete system of differential equations for the 
field amplitude and the collective atomic variables is 
given in Appendix A (|A4|) - (|AlOp . Here, we will write 
the evolution equations of the system in the case where 
the signal field is much weaker than the control field 
|f2p ^ g'^{a^a). Then most of the atoms remain in 
ground state |1) and in the full equation (jA6|) we can 
replace the difference of operators Ni — by the c- 
number N and we can neglect the term proportional to 
(723 • We only keep the coherences <7i2 and aiz- With 
these approximations, one obtain a simplified system of 
the form 



d i d"^ 

'd'z ~ Wd^ ' o,{z,p,t) = ~g ai3{z,p,t), (2.13) 



d_ 

dt 
d_ 
dt 



aisiz, p,t) = gNa{z, p,t) + mi2{z, p,t), (2.14) 



ai2iz,p,t) = -Q*ai3{z,p,t). 



(2.15) 



Due to the assumption that the light pulses are much 
longer than the medium, we have neglected the transient 
regime associated with the time derivative of the signal 

field in Eq. ([203]) . 

We now take the Fourier transform of Eqs. ()2.13p - 
([2l5l) as 

F{z,t;q) = ^ J F{z,p,t)e-''lPd^p, 

F{z,p,t) = ^ / F{z,t;q)e'1P (2.16) 



a(z,t;g)-^a(z,t;q)e-*9'^/(2fc«), 

a^niz,t;q)^ ct™„(z, i; (?) e-*'?'^/(2fc^) . (2.17) 
Then the set of partial differential equations giving the 
evolution of the system reads 



^a(z, t;q) ^ -g aisiz, t; q), 



(2.18) 



—ai3{z,t;q) = gNa{z,t;q) + n&i2{z,t; cf), (2.19) 



— 6-12(2, = ~n*ai3{z,t;q). 



(2.20) 



These equation are similar to the equations obtained in 
Refs. [la, [H, for the plane wave case. However, let 
us underscore that they include the transverse spatial de- 
pendence of the signal field and they allow to treat com- 
pletely the case of a multimode transverse field. More- 
over, we present here a full treatment where the hypoth- 
esis of very large values of ^ gVcN is not made. 



III. ANALYSIS OF THE WRITING PROCESS 

A. Semiclassical solutions of the equations for the 
writing process 

The aim of this section is to derive exact solutions for 
the system of equations (|2.18p - (|2.20p . in order to have a 
detailed information on the efficiency of the writing and 
reading processes in the atomic medium in various condi- 
tions. This will allow optimizing the storage and retrieval 
processes. In order to solve these equations, we will use 
the Laplace transformation in the time domain [l^ [sl] . 
In this case two of the three differential equations are 
transformed into linear algebraic equations, and the so- 
lution of the third differential equation can be written 
in an explicit form. An inverse Laplace transform then 
gives the values of the amplitude a(z, t; q) and coherences 
ai2{z,t;q) and ai3{z,t;q) under arbitrary initial condi- 
tions. The detailed description can be found in Appendix 
A. In the following we will concentrate on the efficiency of 
our memory model, for which the semiclassical solutions 
are sufficient [H, [la] . 

According to Appendix A the semiclassical solutions 
for the writing stage for arbitrary relations between g^cN 
and I rip read 
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a^{t,z;q)^ / dt'a^n{t' ; q) D{t - t' , z), 

JO 

a^{t,z;q) = -gN [ dt' sin \n\{t - t') [ dt"a,r,{t" ;q) D{t' - t" , z), 
Jo Jo 

rt ft' 



ar3{t,z;q)=gN / dt' cos \n\{t - t') / dt" a.,n{t" -q) D{t' ~ t" , z), 
Jo Jo 

where the kernel D{z,t) is expressed via the first-order Bessel function of the first kind Ji in tlic form 



(3.1) 
(3.2) 
(3.3) 



D{z,t) = S(t) - cos\n\t 



2g^Nz 



t 



J I {y/2g'^Nzt 



+ l9^NzJ^ dt' 



l=Sm-t')j^^^2g^Nz{t-V)) 



Vt-t 



(3.4) 



To examine the efficiency of the storage for each spatial spectral field component, we assume that an amplitude of 
the input signal pulse is constant in time, 



In view of numerical calculations, the solutions can then be written in the following form 



a (t,z;q) = ainiq) 



1+ / dt' b{i',z) 



<^12 {i,z;q) = -p Qiniq) 



l-cosf+ [ di' [l - cos{i - i')] D{i',z) 
Jo 



<^i3{i,z;q) =pain{q) 



smt+ I dt' sm{t-i') D{i',z) 





where we define a new kernel D{t, z) as 



D{t,z) = \n\ S{i) + D{i,i) 



and 



D(i,z) = - cost \^Ji(Vzi) + -z di' 
\ t ^ ' 4 Jo 



4=e~'*'ji(yiP) 



Vi - i 



(3.5) 

(3.6) 
(3.7) 
(3.8) 

(3.9) 
(3.10) 



r 



We have defined an effective interaction coefficient p, 
given by 



P ■■ 



gN 



and we have introduced dimensionless values of for time 
t and space z, defined by 



t = \n\ t. 



- 2g'N 
z = ,„, z. 



(3.11) 



The above definitions can be understood in the following 
way. We focus our analysis on a time scale on which the 



spontaneous decay of the upper level is negligible, so the 
effective evolution rate of the levels is determined by the 
Rabi frequency fl. The inverse Rabi frequency Q,^^ is 
then the natural time unit. Moreover, if we replace the 
relaxation constant of the upper level 7 by an effective 
decay rate which is the Rabi frequency in the usual ex- 
pression of the optical depth for an atomic medium of 
length z, we get an effective optical depth 2g'^Nz/\n\. 
This expression is the dimensionless space coordinate z 
defined above. 

The expression for the probe field amplitude (|3.1I) is 
the convolution of the field value at the input of the 
medium with the kernel D{z,t). The expressions for 
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the coherences (|3.2|) - (|3.3|) involve double convolutions. 
Rather than looking for an analytical solution, we have 
solved this system numerically, and in the next sections, 
we analyze the results of this calculation. 



B. Evolution of the signal field and of the atomic 
coherence inside the memory cell 

As stated earlier, we study here the light-matter inter- 
face in conditions where the light pulse is much shorter 
than the atomic excited state lifetime but much longer 
than its propagation time in the atomic medium L/c. 
Thus the propagation time of the pulse wavefronts in- 
side the medium is very short, and we can neglect the 
evolution of the field and the atomic state and also any 
energy exchange between field and atoms over times on 
the order of L/c. In this approximation, we also consider 
the driving field as constant during the pulse duration 
Tw . Note that the latter approximation is not funda- 
mentally necessary for our calculations but it simplifies 
the solutions. 

In view of the above approximations, we can define 
the writing process duration as the interval between time 
t = 0, when the front part of the pulse goes out of the 
medium and time t = Tw when the end of the pulse 
reaches the input surface of the medium. 

Let us normalize the amplitude {i,z,q) of the in- 
put field with respect to its value before it enters the 
medium, and the coherence a^{t,z,q) with respect to 
its maximal possible value at point z — 0, which can be 
shown to be —2pain{q) as f = tt from Eqs. p.6p - (|3.7p . 
The normalized signal field and atomic coherence can be 
written as 

a \t^z) = -— , 

„W (1 ~\ Q'l2 ft; q) 

(^I2yt,z) = -— (3.12) 

2pain(q) 

We first analyze the distribution of the field amplitude 
inside the medium as a function of time. At time t = 
the signal field is distributed homogeneously along the 
medium and is equal to 1. It remains equal to 1 at the 
entrance of the medium, z = till the end of the signal 
pulse at t = Ty^. 

The three upper panels (a, b and c) of Fig. [5] show the 
calculated distribution of field amplitude along z in the 
atomic medium for times t — 0.5, 1 and tt. The amplitude 
of the field shows an oscillatory behaviour as a function of 
z, the effective optical depth at point z. This behaviour 
results from the interplay between the signal field and 
the atoms interacting with the strong driving field. 

Because of the interaction of the atoms with the driv- 
ing and signal fields a small fraction of the atoms ini- 
tially placed in state |1) is transferred to state |2) and an 
atomic coherence ai2(t, z) is generated between states |1) 
and 1 2). The evolution of this coherence is shown in the 



three lower panels of Fig. [2j This long-lived coherence is 
the basic feature of the quantum memory. 

The lower panels in Fig. [5] clearly show the build-up 
of the ground state coherence, which results from the in- 
tegrated interaction of the atomic medium with the field 
since the beginning of the pulse. Our calculation allows 
to follow in detail the build-up of this coherence. Up to 
t = TT (which corresponds to a tt driving pulse), the co- 
herence {i, z) increases dramatically with time up to 
its maximum value of 1 for small values of z and stays 
close to for large values oi z {z > 15). Thus only the 
layers close to the input of the atomic medium are in- 
volved in the storage process, while deeper inside of the 
atomic medium, the coherence remains equal to zero. 

For longer time intervals the behaviour of the atomic 
coherence changes. Figure [3] depicts the evolution of 
a^{i,z) as a function of i and z. At times later than 
t — TT the distribution the coherence over the medium 
changes: decreases for small z and starts to grow for 
larger values of z. It can be seen clearly in Fig. [3] that 
the maximum of shifts from z = to non zero values 
of z. 

Such a behavior is the result of two competing pro- 
cesses: local writing and reading of the signal field. Ac- 
tually, after time t = tt, we can consider that the control 
field starts to "read" the atomic coherence and writes 
it back further in the medium. Quantum information 
is transferred from the input layers to the deeper ones. 
Eventually, the whole medium is involved in the storage 
process. 

If the duration of the field-medium interaction is in- 
creased even further, the maximum of the atomic coher- 
ence shifts deeper into the atomic medium, as shown in 
Fig. 01 for i = TT, 2tt, 3tt and Att. 

Let us note that the behavior of the coherence in a 
thin layer close to the input is similar to the one pre- 
dicted in Ref. [Hj for the case of a single-mode cavity: 
the coherence increases from to 1 between times and 
TT, then it decreases from 1 to in the following tt in- 
terval, etc., that is, the writing and reading processes 
interchange each other with a period tt. However, in our 
case, this is only true for z ~ 0. For other points inside 
the medium, the behaviour is more complex than this 
periodic oscillation. 



C. Optimization of the writing process 

In the previous section, we have analyzed the writing 
process, that is the build-up of the ground state coher- 
ence in the atomic medium. From the calculated dis- 
tribution of fjj^, it can be seen that an efficient writing 
process, with an effective state exchange between signal 
field and atomic coherence can only be obtained for a 
pulse duration larger than tt. However, as explained be- 
low, the TT pulse does not necessarily correspond to the 
optimal duration because it only ensures an effective ex- 
citation of the input layers. In order to quantify the 
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FIG. 3: Distribution of the coherence ai2{i,z) in time and 
space. 



writing efficiency, we will use the intensity of the out- 
going signal field, since any photon from the signal field 
going out of the medium can be considered as a loss. 
Using the method presented in the previous section, we 
obtain the signal field at the output of the medium of 
length L, a^{i,L). The time dependence of this field 
normalized according Eq. p.l2p is shown in Fig. [5] for a 
specific length of the medium. It can be seen that this 
outgoing field, usually called leakage may take non neg- 
ligible values. We will characterize the efficiency of the 



memory by the total losses 

fi^"" \a^(i,L)\^di 
Loss{Tw, L) = ^-^-J ^-^-^ X 100% (3.13) 

For a length L of the medium and a duration of 
the writing pulse, the losses correspond to the integral of 
the signal field intensity over the duration of the writing 
pulse normalized to the full pulse energy. 

Using Eq. ()3.13|) , we can optimize the writing effi- 
ciency by minimizing the losses. This is shown in Fig. [51 
where the upper panels depict the losses as a function of 
the pulse duration for various medium lengths. 

Figure ini shows that for each length, there is an optimal 
pulse duration that minimizes the losses. For a medium 
length of L = 10.3, the minimal losses are obtained for 
TVi/ = 4.2 and are about 5.1 %. The losses decrease with 
longer atomic media, but at the expense of longer writing 
pulses. 

For each of the considered medium lengths, the lower 
panels of Fig. [6] show the distribution of the atomic 
ground state coherence corresponding to the optimal 
writing time Tvv- For short atomic media and pulse 
time durations. Fig. losses are significant due to a 
field-medium interaction that does not last long enough 
for the build-up of an appreciable atomic coherence. On 
the other hand, if the pulse duration is longer than opti- 
mal for a given medium length, signal field leakage takes 
place: the atomic coherence starts to be re-read by the 
control field resulting in a field emission that increases 
the losses. 

The losses Loss(Tw,L) can also be calculated as a 
function of the medium length for a fixed pulse duration. 
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FIG. 5: Field a^{i,L) at the output of the medium during 
the writing process for L = 10.3. 



This is shown in Fig. [71 In this case, the losses decrease 
monotonically when the medium length increases. Thus, 
if the atomic medium can be considered as an unlimited 
resource, any signal pulse (satisfying Eq. (12.11) ) can be 
written with a predetermined efficiency. On the other 
hand, if the length of the medium is limited, then one 
can no longer write an arbitrary pulse and optimization 
is required. 

Finally let us examine how realistic the parameters re- 
quested for optimal writing are. It can be seen in Fig. 
|6}d that the ratio of L/Tw, ensuring minimum losses is 
approximately equal to 2.5. Turning back to dimensional 
variables and separating three main factors, we get 



w 



2g^NL 1 7 
7 jTw l^ip 



TKH^ (3-14) 



where 7 is the spontaneous decay rate of the upper level. 
The first factor is the (real) optical depth, that is typi- 
cally on the order of 1 to 10. The second factor, assuming 
(|2.ip is much larger than 1. Thus, the Rabi frequency 
is determined by the inequalities 7^/|r2p ^ 1. Also one 



can see that these optimal conditions are compatible with 
g^/Nc ~ |ri|. The conditions necessary to reach a good 
efficiency are expected to be quite feasible. 
IV. ANALYSIS OF THE READ-OUT PROCESS 

A. Semiclassical solution of the equations for the 
read-out process 

The quantum information written into the atomic 
medium as described in the previous section can be stored 
during a time which is limited by the the lifetime of the 
ground state atomic coherence. For times shorter than 
this lifetime, the quantum information can be recovered 
by means of the read-out process. The latter is carried 
out using a driving pulse at the same frequency as the 
one in the writing stage. We will consider two possible 
geometries for the read-out pulse propagation. The first 
one is the forward retrieval for which the read-out driv- 
ing pulse propagates in the same direction as the writ- 
ing pulse. The second one is the backward retrieval for 
which the read-out pulse has the opposite direction. Like 
for other memory schemes [3 l39l442j , we will see that 
backward retrieval provides better read-out than forward 
retrieval for a single mode while it is not necessarily true 
for a multimode signal. 

The solutions of equations (|2.18p - (|2.20p were obtained 
in the previous section for specific initial conditions char- 
acteristic of the writing process. For read-out it is as- 
sumed that the signal field at the input of the medium 
a^{t,z = 0;q) and the coherence at the initial time 
a^^{t = 0, z;q) are in the vacuum state, i.e., equal to zero 
in the semiclassical approach. The ground state coher- 
ence a'^2it — ^1 ZjQ) is assumed to coincide with its value 
at the end of writing stage. First in the case of forward 
read-out, the solutions describing the reading process are 
given by (see Appendix A) 
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FIG. 6: Writing process : relative losses in the field intensity (in percent of the input field intensity) at the output of the 
medium as a function of Tw for (a) L — 5; (b) L — 10.3; (c) L = 15 and corresponding coherence distribution for optimal 
conditions. 
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FIG. 7: Writing process : relative losses in the field intensity 
(in percent of the input field intensity) at the output of the 
medium as a function of L for Tw ~ ti". 



a^{t,z;q)^-g [ dz f dt' sm\n\t' aY2{Tw , z' -q) D{t - t' , z - z'), (4.1) 
Jo Jo 

(^i^it, z;q)^ sin \VL\t (J^{Tw, z;q)~- 

-g^N [ dt' cos\Q.\(t-t') [ dt"sm\n\t" [ dz'a^{Tw, z' ; q)D{t' - t" , z - z'), (4.2) 
Jo Jo Jo 

a^^it, z- q) = COS \n\t ai^(T, z- q) + 

+g^N [ dt'sin\n\{t-t') [ dt" sin \n\t" [ dz'a^{Tw, z';q)D{t' ~ t" , z ~ z'), (4.3) 
Jo Jo Jo 

I 



where the ground state coherence at the beginning of 
the process cr]^(Tvi/, z', q) contains the stored signal field. 
This value coincides with equation p.2p at t — Tw- 



Solutions for the backward read-out are similar except 
for the spatial argument in the coherence a^{Tw,z';q) 
where z' has to be changed to L — z'. Substituting the 
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coherence ai2 from Eq. (|3.2p to Eq. (|4.ip and taking into 
account the diffraction factors (|2.17l) one can derive the 
signal field at the output of the memory cell for the for- 
ward read-out of the form 



1 







sint+ / dt'sini' D{i-i\L-z') 
Jo 

ffw 

1 - cos fw+ / di" 
Jo 



X 



(4.4) 



l-cos{Tw-i") D{t'\i') 
and correspondingly for the backward read-out: 



focus on restoring the pulse shape, but rather, on the 
read-out of an image corresponding to an ensemble of 
the transverse modes of the signal field, which has been 
stored in the atomic medium as a quantum hologram. To 
do this, the photon number in each of the reconstructed 
transverse modes must be the same as the photon number 
in the corresponding modes of the original field. To esti- 
mate the storage efficiency, we introduce the parameter 
Eff{q), defined as the ratio of the integrated intensity of 
the recovered pulse in mode q to the integrated intensity 
of the same transverse mode in the input pulse : 



Effiq) = 



\a''{t,z,q)\^dt 



/o \ain{i,q)\'^di 



X 100% 



(4.6) 



''hack 



(i,L;q} = -a,^{q) / dS' e-'"" Z"^' 



sin t -f / di' sin i' Dii-i', z') 
Jo 



l-cosT, 



w 



Tw 



dt'- 



(4.5) 

l-cos{fw -i")] b{i",z') 



These equations are written using the dimensionless time 
and space coordinates except for the exponential diffrac- 
tion factor where the regular z coordinate is kept. One 
can see that the situation is very different for forward and 
backward retrieval. For forward read-out the diffraction 
introduces a common factor which can easily be compen- 
sated with an appropriate lens. Thus, in this case, al- 
though diffraction takes place separately in writing and 
reading, the effects of diffraction due to the writing pro- 
cess can be perfectly compensated under the read-out. 

A dramatically different situation takes place for the 
backward read-out. Because the diffraction factor is un- 
der the integral, the diffraction are able to modify se- 
riously the initial transverse distribution. In contrast 
to the first case, the diffraction processes in the writ- 
ing and reading channels do not compensate each other 
but on the contrary add their effects. From Eqs. (14. 4^ - 
(|4.5[) . it can be seen that under the geometrical condi- 
tions q^L/ks ^ 1, diffraction is not detrimental (see also 
below). Here we shall consider only this case and we 
show in Fig. 8 the time dependence of the intensity of 
the read-out signal. 



B. Efficiency of the memory process 

Fig. |S] shows the intensity of the re-emitted normalized 
signal field a^{i,L) — a^{i,L,q)/ain{q} at the output of 
the medium as a function of the reading time for forward 
and backward retrieval for an optical depth L = 10.3 and 
a writing time duration Tw = 4.2. It can be seen that 
the shape of the restored pulse is quite different from the 
shape of the input signal pulse. However, we will not 



where Tr is the reading time duration. It can be seen 
from Fig. |5] that for an efficient recovery, for forward 
retrieval, the reading time duration should substantially 
exceed the writing time duration. For example, in Fig. 

the retrieval efficiency after a reading time equal to 
the writing time Tw is only 36%. For backward retrieval, 
the retrieval is much more efficient. In Fig. [5]d the effi- 
ciency reaches 80% for a reading time Tr — Tw- 

Losses come from both the writing and the reading 
processes. With parameters L = 10.3 and Tw = 4.2 Fig. 
[S] shows that the losses in the writing process are 5.1 %. 
Then, the best efficiency of the whole process, including 
writing and read-out cannot exceed Eff = 94.9%. 

For the reading process, a numerical calculation shows 
that if the reading time is 10 times the writing time 
Tr = lOTw, the overall efficiency is Eff = 77% for 
forward retrieval. We could choose an even longer read- 
ing duration, but we must remain in the approximation 
that Tr ^ 7~^. The backward retrieval is faster and the 
same calculation yields an overall efficiency of 83.8% for 
Tr — 3Tw- 

In the same way as for writing, we observe oscillations 
in the read-out field in forward direction (Fig. [8^). This 
indicates that the reading process is accompanied by a 
rewriting process in the medium, before emission. In the 
case of backward retrieval, the situation is much more 
favorable and the curve does not show any oscillation 
(Fig. Eb). 

The read-out as well as the writing is coupled with a 
displacement of the coherence along the medium. From 
the results obtained in the previous section it can be seen 
that for not too long writing pulses, the storage in the 
atomic coherence is mainly localized in the input layers 
of the medium. Then, if the medium is long enough, the 
forward retrieval involves reading the signal and writing 
it again on successive layers in the medium, which even- 
tually "pushes" the read-out signal out of the medium. 
On the contrary, in the backward retrieval there is less 
" pushing" , since a large coherence is concentrated at the 
exit face. 
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FIG. 8: Reading process : field intensity \a^{i,L)\^ at the output of the medium for L = 10.3 and TV = 4.2 for (a) forward 
and (b) backward propagating retrieval. 



C. Restoring the transverse profile 

Let us now examine in more details the information 
about the transverse profile of the field in the recon- 
structed pulse. Usually, the quantum memory of the 
pulse implies its longitudinal profile storage. A lot of 
works devoted to this problem (for example, [? ].(l6|.[l8|) 
have focused on the determination of eigenfunctions of 
the system, we will not deal with this issue. Our goal 
is to restore the transverse profile of the signal field. To 
achieve a quantum memory we need to ensure high ef- 
ficiency of reconstruction for each transverse mode. In 
the previous section, we showed that each of the trans- 
verse modes can be recovered with equal efficiency, with 
a value that depends on the choice of experimental pa- 
rameters and that can be high enough. However for the 
quantum storage of the transverse profile we should have 
not only a high efficiency for each component but also an 
insensibility to the diffraction phenomena. 

In the previous section we showed that the reading 
conditions are different depending on whether reading is 
a direct or reverse process. However, this apparent ad- 
vantage of backward reading is associated with a limita- 
tion due to diffraction phenomena. The essential differ- 
ence between the forward and backward reading is that 
in the first case, the diffraction distortions that occur 
during storage are almost entirely compensated under 
subsequent reading. What is left uncompensated can be 
easily cleaned, for example, by setting additional lens at 
the output of the cell. 

In the second case, the diffraction distortions that oc- 
cur sequentially during recording and reading, add to 
each other, and they have to be minimized independently. 
It can be seen that this leads to a significant limitation for 
an important memory parameter, which is the number of 
modes or the grain of the image that can be efficiently 
stored in the memory cell. In principle, this number of 
modes is determined by the ratio between the transverse 
size of the memory cell S and the grain size d of the im- 
age at the entrance to the cell and is equal to = S/(P. 
Minimizing the diffraction means that at the output of 



the memory cell of thickness L the grain size D = LX/d 
should coincide with the input grain size D ^ d. This 
means that the number of stored modes N can not exceed 
the Fresnel number Fjy — S/{XL). 

So our numerical analysis of backward reading that 
does not take diffraction into account can be considered 
as justified, provided that the number of stored modes is 
below the Fresnel number N < F^. 

Similarly, we can estimate the number N for forward 
reading. In this already mentioned, we can gen- 

erally ignore the diffraction image distortion associated 
with the intersection of various grains at the cell output. 
However, we should not allow excessive beam divergence 
and, as a result of this divergence, light losses due to 
leakage of the field through the side walls of the cell. To 
avoid these losses, the output grain size should not ex- 
ceed the transverse size of the memory cell D < y/S. It 
is easy to see that the number of stored modes should 
not exceed the square of the Fresnel number N < F^, 
in agreement with Ref.fsl]. We can thus predict a large 
storage capacity for this A-based memory scheme. 



V. CONCLUSION 

We have presented a full calculation of the writing and 
read-out processes of a short pulse in an atomic medium. 
Our approach relies on a full treatment of light matter 
interaction as the writing or reading pulses propagate 
through the atomic medium. This allows to obtain de- 
tailed predictions on the behavior of the atomic ground 
state coherence, which stores the quantum information 
and on the re-emitted field. As a result, precise values 
of the storage efficiency and of the memory capacity for 
multimode storage are given, together with a detailed op- 
timization procedure for the memory operation. We have 
also shown that if backward retrieval provides a better 
efficiency for a single mode, in the case of spatial multi- 
mode storage it tends to limit the number of stored modes 
to a lower value than the forward retrieval. With the ex- 
perimental developments of short pulse storage and fast 
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memory registers, our model provides a reliable method 
to optimize the experimental conditions. 
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Appendix A: Main equations and general solutions 

The light matter interaction Hamiltonian for our prob- 
lem can be written as 



V 



dz d 



■p[ig (a(z,p,t)e*^-^ ^^h^i{z,p,t)- 



tn{t)e'''d'^-'^hs2{z,p,t)- 
^n*{t)e-'^'^' + '^h,:,{zJ,t) 



At- 



(Al) 



where fc^, kd are wave vectors of the signal, driving waves, 
and the value A determines the two-photon resonance 



A = CJ^ - Wis = - UJ23- 

We perform the substitutions 



(A2) 



- , JkgZ — iAt- 
(713 -> e ^13, 



- , JkdZ — iAt ' 

CT23 e " CT23, 



0'12 



'i{kd - ks)z 



0'12. 



(A3) 



Then the Heisenberg evolution equation for the atom field 
system can be written as 

d d ^ 

dt dz 2ks 



d_ 
dt 
d_ 

di 
d_ 

di 
d_ 

di 
d_ 

di 



0-13 



0'12 



= -cgai3, (A4) 
iA(7i3 + flai2 + gaiNi - N3), (A5) 
n*ai3~gda32, (A6) 



6-32 = iAo-32 - n*{N3 - N2) + ffatai2, (A7) 



N2 



dt 



-gaa^i - ga'a^, 



~f7(732 — il*(723, 



(A8) 
(A9) 
(AlO) 



where the transverse Laplace operator reads 

92 92 



A, = 



dx^ dy'^ 



(All) 



Here we do not take into account a spontaneous emis- 
sion with rate 7 on the transition |3) — > |1). This is 
coupled with a requirement that in our model the field 
pulses are very short such that the spontaneous emission 
has no time to introduce something to the atomic state. 

By ignoring in Eqs. (jA4p - (jA10[) the spontaneous relax- 
ation we nevertheless survive the frequency detuning A. 
The condition A ^ 7 usually means that we want to 
treat the Raman process. However for the short pulses 
this condition do not ensure Raman process, here Raman 
interaction is realized only in the case of stronger inequal- 
ity: A ^ T~^. In this article we consider the opposite 
condition A <C T^^, when we are under obligatory have 
to neglect by A in Eqs. (KE\i and (K7\ . 

From Eqs. (|A4[) - (|A10p one can obtain that the operator 

S{t) = f d\ (h{r, t)lc + N2{r, t) + N^if, t)) (A12) 



survives in time. Here n{r,t) = a) {r^t)d{r,t) is the pho- 
ton number operator per second per unit area. Physically 
this result seems quite natural. 

We assume the following conditions. All the N 
atoms are initially in the state |1) and we assume that 
the signal field is much weaker than the driving field 
I rip 3> g"^ {a)a). Then most of the atoms remain in 
the ground state |1) and in Eq. (jASP in approximation of 
the quasi-homogeneous medium we have a right to make 
a change (iVi ~ N^) — > N and also we can neglect the 
term proportional to (T32 in Eq. (jA6p . Thus we obtained 
a reduced closed system of equations for interesting val- 
ues 



^ d d i . , ^ , _ , 

c9^ + 9;-2fc:^^'"^"'^'') 

'9 cri3{z,p,t), 



(A13) 



d_ 
di 
d_ 

di 



^isiz, p,t) ^ gNd{z,p,t) + nai2{z,p,t), (A14) 



^i2iz,p,t) = ~n*ai3{z,p,t). 



(A15) 



We neglect the time delay linked to the pulse propaga- 
tion in the atomic medium. This means if we have long 
enough pulses, such that L/c^T (L is thickness of the 
medium and T is the pulse duration) , we can neglect the 
time interval between the time at which the front part of 
the pulse enters the medium and the time at which the 
front part leaves it. Formally this means we can neglect 
the time derivative in Eq. (jA13p . 

Let us take the Fourier transform and next the Laplace 
transform of the equations according to the relations 



F{z,p,t)e-"1P d^p, (A16) 



Fsiz;q) 



dtF{z, t;q} e 



(A17) 
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In the Fourier domain the equations read 
d 



dz 
d_ 
dt 
d_ 
dt 



a{z,t]q) = -g ai3{z,t;q), 



(A18) 



d-i3{z,t;ci) = gNa{z,t;q) + ilai2{z,t; cf), (A19) 



(A20) 



where we have made the changes 

a{z,t;q) ^ a{z,t;q}e-^1^^/(^''-), 

amn{z, t; q) ^ ct„„ (z, t; q) e-'l''^/ (^fc,) (^21) 

In the Laplace domain the equations are written in the 
form 



das{z; q) 
dz 



(A22) 



-~^i3{z, 0;q) + S(5-i3,s(2; q) gNa^iz; q) + 
f]<7i2,,(z;<7), (A23) 
-<Ti2(2;,0;g) + sai2,siz;q) = 

-n*ai3Az;q)- (A24) 

Eliminating the coherences (Ti3,s and 1T12.S we obtain a 
differential equation for the Laplace field amplitude as : 



ddsjz; q) 
dz 



= -jsds{z;q) - gAs{z;q), (A25) 



where 



7s = g'N 



Asiz;q) 



s2 + |r!|2' 
1 



(A26) 



[fi(Ti2(0, z;q} + S(Ti3(0, z; q)] 



The solution of equation (|A25|) reads 



ds{z:q) = ds{0;q)e 



-IsZ _ 



'AJz':n)p-^s{z-z')_ 



g I dz'As{z';q)e 
JQ 



(A27) 



Taking the inverse Laplace transform we obtain the field 
amplitude expressed as a function of the initial conditions 



a{t,z;q)^ / dt'cimit - t' iq) D{t' , z) 



g dt' dz'A{t-t',z-z';q)D{t\z'). (A28) 
Jo Jo 



The kernel D{t, z) is expressed from the Bessel's func- 
tions and the function A{t, z; q) reads 



t) = 6{t) - cos \Q.\t \ — Ji y/2g^Nzt 



2g^Nz 



t 



^g^NzJ^dt' -^e-*l"l^'ji (TV^) -^=i=e*l"l(* " J, (^^2g^Nz{t - t')) 



(A29) 



A(t, z;q) = cos \Q\t 0-13(0, z;q) + sin \Q\t 0-12(0, z;q). 
From Eqs. (jA19p and (|A20p we also get explicit solutions that read 

ai3{t, z;q) ^ gN [ dt' cos \n\(t ~ t') a{t' , z;q) + A{t, z; q), 



ai2{t,z;q) = e 



-gN / dt' am\n\{t~t') a{t',z;q) + B{z,t;ci) 
Jo 



where 



B{z, t\q) = - sin \Vt\t 0-13(0, z;q) + cos \Vt\t 0-12(0, z;q). 



(A30) 

(A31) 
(A32) 

(A33) 



Thus we have obtained the solutions of the system of the m ain equations under the arbitrary initial conditions. 
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